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ABSTRACT 

GRB 021206 is one of the brightest GRBs ever observed. Its prompt emission, 
as measured by RHESSI, shows an unexpected spectral feature. The spectrum 
has a peak energy of about 700 keV and can be described by a Band function 
up to 4.5 MeV. Above 4.5 MeV, the spectrum hardens again, so that the Band 
function fails to fit the whole RHESSI energy range up to 17 MeV. Nor does 
the sum of a blackbody function plus a power law, even though such a function 
can describe a spectral hardening. The cannonball model on the other hand 
predicts such a hardening, and we found that it fits the spectrum of GRB 021206 
perfectly. We also analysed other strong GRBs observed by RHESSI, namely 
GRBs 020715, 021008, 030329, 030406, 030519B, 031027, 031111. We found that 
all their spectra can be fit by the cannonball model as well as by a Band function. 

Subject headings: gamma rays: bursts — gamma rays: observations — tech- 
niques: spectroscopic 



INTRODUCTION 



The exact mechanism which produces 7-ray bursts (GRBs) has not yet been defini- 
tively established. Their prompt 7-ray spectra can be used to distinguish between different 
models. Several mathematical functions have been used for parametrizing the prompt 7-ray 
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emission. Most commonly used is the empirical Band function (IBand et al.lll993l ). which is 
not motivated by a physical model. 

There have been attempts to distinguish b e tween spectra l mode ls analysing the low 



energy part of t he sp ectrum. iGhirlanda et al.l (120031 ) . iRyde I (12004 ). and more recently 



Ghirlanda et al.l (120071) searched for b lackbody components in GRB spectra with varying 
degrees of success. IPreece et al.l ( 120021 ). using BATSE GRB spectra, tested the synchrotron 
shock model and conclude that it "does not account for the observed spectra during the 
GRB phase". 

Spectral studies above the peak energy are rare, one reason bein g the poor data quality 
because of lack of statistics. Combining BATSE and EGRET spectra, I Gonzalez et al.l (120031 ) 
report a high energy component for GRB 941017. They find a photon index of about 1.0 at 
energies above 5 MeV. 



In this paper we report a high energy component in GRB 021206 (IHurley et al.ll2002c 



2003br). observed with the Reuven Ramaty High Energy Solar Spectroscopic Imager RHESSI 
( iLin et al.ll2002l ). Having a peak energy of about 700 keV, the spectrum of this burst can be 
described by a Band function from 70 keV up to 4.5 MeV, with a high energy photon index 
(3 ~ 3.2. Above 4.5 MeV, the spectrum hardens again, and can be described with a photon 
index f3' ^ 2.2. This significant hardening around 4.5 MeV can not be described with a 
Band function. But it seems to differ from the spectral hardening in GRB 941017 as well. 

T here is one model that fits the entire RHESSI spectrurn of GRB 021206: the cannonball 
model (jPar fc de Ruiulall2004j : iDado. Par fc de Ruiulall2002l . l2003al ). The cannonball model 
predicts a spectral hardening at several times the peak energy with a high energy photon 
index reaching (3 ^ 2.1. 

The question immediately arises whether the cannonball model can improve our descrip- 
tion of other GRB spectra. The difference between the Band function and the cannonball 
model arises only at the high energy part of the spectrum, where data usually suffer from 
low statistics. Therefore, we choose the strongest GRBs registered by RHESSI in the years 
2002 to 2004. We find that they all can be fit by the cannonball model as well as by the 
Band function. 

The outline of the paper is the following: We first present shortly the instrument, the 
GRB selection, the spectrum extraction, and the fit method (§2]). In the next section (^, 
many spectral functions are given. In ^ the fit results for GRB 020715, GRB 021008, 
GRB 021206, GRB 030329, GRB 030406, GRB 030519B, GRB 031027, and GRB 031111 
are presented. The fits are discussed and, if possible, compared to other measurements. The 
more general discussion, including an outlook, follows in ^ We end with a short summary 
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2. INSTRUMENT AND METHOD 
2.1. Instrument 

RHESSI is a NASA Small Explorer mission designed to study solar flares in hard X- 



rays and 7-rays (ILin et al.ll2002l ). It consists of two main parts : an imaging system and the 



spectrometer with nine germanium detectors (jSmith et al 



2002). The satellite always points 



towards the Sun and rotates about its axis at 15rpm. The Ge detectors are arranged in a 
plane perpendicular to this axis. 

The shape of the detectors is cylindrical with a height of ~ 8.5 cm and diameter of 
7.1cm, and they are segmented into a thin front (~ 1.5 cm) and a thick rear segment 
7 cm). Since the shielding of the rear segments is minimal, photons with more than 
about 25 keV can enter from the side. Above about 50-80 keV, photons from any direction 
can be observed. Each detected photon is time- and energy-tagged from 3 keV to 2.8 MeV 
(front segments) or from 20 keV to 17 MeV (rear segments). The energy resolution is ~ 3 keV 
at 1 MeV, and the time resolution is 1 fis. 

The effective area for GRB detection depends on the incident photon energy E and 
the angle between the GRB direction and the RHESSI axis, the incoming angle 6. Over a 
wide range of E and 6, the effective area is around 150 cm^. The sensitivity drops rapidly at 
energies below ^ 50 keV. 



2.2. GRB selection 

For this study, we need well observed GRB spectra. We chose GRBs from the years 
2002 to 2004, because radiation damage starts to play a role in 2005. The selected GRBs 
have to be localized by other observations of the same GRB (RHESSI can not measure the 
incoming angle), because 6 enters into the simulation of the response function. A further 
requirement was the availability of good background data. And flnally, the data storing 
mode ('rear decimation' for onboard memory saving) is not allowed to change during the 
entire GRB and background time interval. Of all the GRBs meeting these criteria, we chose 
eight with the best signal-to-background ratio, listed in Table [1] along with their incoming 
angle 6 and the time intervals used. The lightcurves of these bursts are shown in Figs. [1] and 

m 
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2.3. Preparation and fit of RHESSI spectra 

The method of analysing RHESSI GRB spectra will be described in detail in a separate 
article (E. Bellm, C. Wigger et al., in preparation). 

For each GRB and detector segment, the total spectrum (GRB plus background) during 
the burst was extracted, as well as the background spectra during two time intervals before 
and after the burst. The background was linearly (sometimes quadratically) interpolated 
and subtracted. The exact time intervals are listed in Table [H Then we added all rear and 
all front segments, except for detector #2 which is slightly damaged and has a bad energy 
resolution. Since all GRBs in this study are strong, the observational errors are dominated 
by the statistical error of the GRB counts, not of the background. 



We simulate RHESSI using GEANT3 flCERNlll993[ ). Knowing the direction of the GRB 



from other instruments, we simulate RHESSI's response to photons coming from angle 9. The 
energy of the incoming photons is simulated as a power law spectrum (i.e. dN/dE oc E^'^"^"^) 
with typically ''jsim = 2 . This power law simulation is only a rough approximation and is not 
not intended to represent the intrinsic GRB spectrum, but instead provides simulated data 
representing RHESSI's conversion of photons to counts. The true GRB source spectrum 
is determined via weight factors for the resulting simulated count spectrum, as described 
below. The upper energy limit of the simulated photon spectrum is typically 30 MeV, in 
the case of GRB 021206 even 40 MeV or 50 MeV. This is important, because an incoming 
photon of e.g. 25 MeV may well make a signal of 15 MeV. Rotation angles are generated 
uniformly, i.e. we compute a RHESSI-spin averaged response function. Since the detector 
arrangement shows an approximate 120 degree symmetry, the averaging gives good results 
as long as the analysed time interval is at least one third of the rotation period {Trot = 4 s). 
This was also confirmed by tests. 

The output of the simulation is an event list, or rather a hit list, consisting of all signals 
registered in the Ge detectors. The simulated hit list, having Ng entries indexed by the 
letter /, contains the deposited energy (Ej^^^) as well as the initial photon energy {El"^). The 
measured hit list contains only the observed energy. 

For spectral fitting, the observed energy histogram is compared with a histogram ac- 
cumulated from the simulated hit list. More precisely: The measured histogram can be 
represented by a fc-element vector M with errors Bm-, and energy boundaries i?Q, E\, E2, 
E^. We normalize the histogram M to the total number of counts in the fit range: 
m = M/Cm and am = ^m/Cm, where Cm = -^i the sum goes only over the bins 
included in the fit, i.e. I = {i \ bin i is included in the fit}. The 'theoretical' histogram S 
is accumulated from the simulated hit list. Each entry is weighted with a factor in order to 
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scale from the simulated power law to the probability density which would be expected, had 
we actually simulated the GRB source spectrum dN/dE. The jth bin contains therefore the 
weighted sum of all simulated hits with Ef'^^ belonging to that bin, i.e.: 



S, = J2w, (1) 



leL 



where L = {I \ E'!_^ < Ef^^ < E°} and 



The first factor in eq. [2] accounts for the spectrum assumed in the simulation and the energy 
Epiy is an arbitrary normalisation. The second factor accounts for the spectrum of the 
incoming GRB photons. Possible parametrisations of dN/dE are given below in ^ If the 
GRB spectrum had the same shape as the simulated one, i.e. if dN/dE = [E / Epi^)''^"^"' , 
the weights would all be 1 . This me thod of using weigh t facto rs when filling a histogram is 
common in particle physics, see e.g. iBarlow fc BeestonI Jl993h. The statistical error of the 



'theoretical' histogram S is cr|^ = '^ifziwf (§6 of iBarlow fc BeestonI Il993l ). As in the case 
of the measured histogram, the histogram S is normalised: s = S/Cs and as = as/Cs with 

— * 

The parameters of the histogram S are varied until the minimum of 



2 (jUi-fSiY 
2^ rr2 Mrr2 



is found. The factor / accounts for the normalisation between measured and simulated 
histogram. It is expected to be almost 1, but should be treated as a free fit parameter. For 
each fit iteration, the histogram 5* is recalculated with different weights (eq. [2]). 

Since the simulated hit list contains many more photons than the measured spectrum, 
we used the approximation a"^. + a^. ~ cx^. while fitting. But it was always checked that 
the statistical error from the measurement is dominant. 

It is possible to create a response matrix from our simulations and perform spectral 
fits via forward-fitting, as in XSPEC. In any case, our weighted histogram method gives 
equivalent fits to response matrices which are simulated directly (by EB; see E. Bellm, C. 
Wigger et al., in preparation). 
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2.4. Systematic effects 

At low energies, a small deviation of our RHESSI mass model from the true amount of 
material can make a considerable difference in the number of observed photons. For 9 ^ 90 
degrees, this should be a small problem because the lateral shielding is thin. But for 9 from 
10 to 50 degrees this is an issue, and less prominently also from 130 to 160 degrees. 

Simulation quality also gets better with higher energy. This is fortunate for the current 
analysis which relies on high-energy properties of GRB spectra. 



3. SPECTRAL MODELS 

Let dN/dE be the number of GRB photons per energy bin. The peak energy -Epeafc is 
defined as the energy for which E-Fe = dN/dE E'^ is maximal. The spectrum in the EEe 
representation has at least one maximum because the total emitted energy must be finite: 
Jq°° dN/ dE E dE < oo. Many instruments can see such a maximum Epeak within their energy 
range. 

Different mathematical functions, sometimes called models, can describe such a shape. 
A collection is presented in this section. 

The simplest spectral function is a power law (PL): 

where Ep^y is a normalization energy, e.g. Ep^^ — 100 keV. The PL has no peak energy. It 
rarely fits a GRB spectrum over the entire observed energy range, but it is often useful for 
a limited energy band. Indeed, every spectrum can be fit by several joined PLs. 

One simple way to account for a spectral softening and a peak energy is the cut off 
power law (CPL): 



dE \E ^ 
If a < 2.0 then Ep^ak ^ Eo{2-a) . 

Another way to account for a spectral break is the broken power law (BPL) consisting 
of two joined PLs 



dN _ J A{f) iiE<Eb 



E^^P ..... (6) 



dE {fY iiE> E,. 

If q; < 2.0 and /3 > 2.0, then Ep^ak — Ef, . This function is not continuously differentiable 
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A s mooth transition b etween the two power laws is reahzed by the empirical Band 
function (IBand et al.lll993l ). This is a smooth composition of a CPL for low energies and a 
PL for high energies: 



dN 
dE 



break 

if E > Ebreak 



(7) 



where 



and 



Again, 



Ebreak = Eo{P — a) 



B = A 



En 



E, 



a 



-(/3-a) 



Epiv is a normalization energy, e.g. Epi^^ = 100 keV. If a < 2.0 and /? > 2.0, then 



a 



If p 



oo or if Ehreak hes at the upper limit of the obs erved energy 



range , then the Band function turns into a CPL. As already pointed out by iPreece et al. 



( 120001 ). section 3.3.1., the low energy photon index, the curvature of the spectrum, and its 
peak energy are represented with only two parameters, a and 

Sometimes a blackbody spectrum plus power law is us ed for spectral fitting, see e.g. 



Rvdel t004 \: iGhirlanda et all (120071 ): iMcBreen et al.l (120061 ). We will call this the BBPL 



{E/{kT)f 



- A 

dE ~ exp{E/{kT)) 



+ Ab 



E.. 



piv 



E 



We choose Epi^ = 3.92 kT, the peak energy of the blackbody component. The BBPL function 
can fit a spectral hardening at high energies. 

When fitting with the BBPL model, it is often found that the PL c ompone i it doe s not 
fit simultaneously at low and at high energies. This is also mentioned by iRyde I (120041 ). We 
therefore invented a blackbody plus modified power law (BBmPL): 



dN {E/{kT)f 
dE exp{E/{kT)) - 1 

A6(l-e-^/^o) 

\ E 



+ 



(9) 



^The smoothly broken power law model (SBPL, see IPreece et al.l (|2000l ): iKaneko et al.l (j2006l )) would 
account for this problem with an additional parameter, but we do not use it here, because it can not fit a 
spectral hardening at high energies. 
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We choose again Eq = 3.92 kT. The modification of the power law component was bor- 
rowed from the cannonball model (see next). The BBmPL function can describe a spectral 
hardening at high energies. 



The cannonball model CBM flDar fc de Ruiulal l2004l : IPado. Par fc de Ruiulal 12002 



2003al ) makes a prediction for the spectral shape of the prompt GRB emission. It consists 
of a GPL and a modified power law: 



(10) 



according to lDar fc de Rujulal (l2004j ). eq. 47, or lDado. Par fc de Rujulal (120041 ). eq. 13. The 
theoretically expected values are a ~ 1.0 and jS ~ 2.1. The CBM function eq. [TOl applies 
strictly speaking only to the spectrum caused by a single cannonball, i.e. for every single 
peak of a GRB. 

It is often observed that the peak energy Ep^ak is a more stable fit parameter than 
the parameter Eq in the Band function (eq. [7]) or in the CPL (eq. [5]). Therefore, we use 
Ep = Eq{2 — a) as a fit parameter. Similarly, we use Tp = T {2 — a) as a fit parameter in 
the case of CBM (eq.OTO]). 

A word about fitting CBM: For the high energy part, it has two parameters, whereas 
the Band function has only one. Already when fitting the Band function, it is often observed 
that the high energy power law index is poorly constrained, because the high energy data 
tend to have large statistical errors. This is even worse for CBM with two high energy 
parameters. It often helps to freeze the parameter (3 at its theoretical value of /5 = 2.1 in 
order to make the fit converge. 



4. FIT RESULTS AND FIT DISCUSSIONS 

The spectral models used and the of the fits are listed in Table [2] for all eight GRBs. 
For CBM and Band function, the fitted parameters are listed in Tables |3] and HI respectively. 
Throughout this article, all errors are symmetric la errors if not stated otherwise. 

The measured spectra together with the CBM and Band fits are shown in Figs. [3] to 
Uni Note that we plot energy^ ■ counts/keV versus energy. The difference to a deconvolved 
EFe distribution is discernible e.g. in the drop of counts towards lower energies in our repre- 
sentation. The statistical scatter from the limited number of simulated events is sometimes 
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visible as a little roughness of the simulated spectrao 

From the fit parameters obtained for the CBM and the Band function, we calculate the 
fluences for various energy intervals. They are listed in Table |3 The error of the fiuence 
is dominated by systematics, e.g. because we do not know the exact active volume of the 
single detector segments. We estimate the systematic error to be of order 5%, whereas the 
statistical error is of order 1%. Note also, that the two fluences obtained by fitting CBM 
and Band function are nearly equal. 



4.1. GRB 020715 

The lightcurve is shown in Fig. [1] (top) and the spectrum is shown in Fig. [3l Two bins 
from 290 keV to 310 keV (a background line) and one bin from 500 keV to 525 keV (the 
511 keV line) were omitted in the fit because they would dominate . The best fit is a Band 
function, but the CBM also fits well. 



4.2. GRB 021008 

Coming from a direction about 50 degrees from the Sun, this GRB deposited photons 
not only in rear detectors but also in the front detectors, as can be seen from the lightcurve 
in Fig. [H The front spectrum is shown in Fig. H] and the rear in Fig. O 



Fit We had difficulties to fit front and rear spectra consistently below 300 keV. We there- 
fore chose 300 keV as the lower energy bound. Band function and CBM give the best fits. 

We fitted front and rear segments separately, as well as jointly. The results of the joint 
fit are shown in Figs. |l]and[5] . For the joint CBM fit we find the 90% confidence level (CL) 
errors: 

Tp = 641+^^keV 

a = 1.523 1°^^ (11) 
f3 = 2.1 



^ In the case of e.g. GRB 021206, rear (see Fig. [7]), the mean measured error between 4 and 12 MeV is 
0.65 • 10^ counts-keV, whereas the mean scatter of the simulated histogram is 0.25 • 10^ counts-keV. For the 
other GRBs with less observed photons and therefore larger measurement errors, the statistical error of the 
simulation is even more negligible. 
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h = 0.0198 IHJ^? 

The total is 71.2 for tidof = 74. The Band function fits marginally better, = 70.0 for 
TT-DoF = 74 and its parameters are (90% CL errors): 

Ep = 677 til keV 
a = 1.493 (12) 
f3 = 3.73 



Discussion We do not well understand the spectrum below 300 keV. Both fits, the CBM 
and the Band function, overestimate the counts below 300 keV. This can be a hint that the 
GRB spectrum hardens below 300 keV. We find functions that fit the front and the rear 
data from 40 keV to 400 keV individually, but they do not agree. 

One possible explanation is the GRB incoming angle of about 50 degrees at which the 
GRB photons pass through a certain amount of material before reaching the detectors. Our 
GEANT simulation tries to take that into account, but it is probably not perfect, and maybe 
the averaging over all rotation angles is a bad assumption for this short GRB pulse. 

Another difficulty for this GRB is its background. For the single rear segments, the 
background at low energies (below ~ 120 keV) strongly depends on the rotation angle of 
RHESSI. We did our best to take this into account, but maybe did not succeed completely. 



4.3. GRB 021206 



GRB 021206 is famous for its cla imed polarization (ICo 



Durn 



ever turned out to be an artefact (see 



GRB was also studied by iBoggs et al 



Rutledge fc Foxll2004 



fc Boggsl[2003[) whic h how- 



Wigger et al.ll2004l . l2005f ). This 



(120041 ) to probe quantum gravity. 



This GRB is only 18 degrees from the Sun, exposing mainly the front segments of 
RHESSI's detectors. Its lightcurve is shown in Fig. [H Fig. E] shows the energy spectrum in 
the front segments. Fig. [7] in the rear segments. 



Fit The front spectrum can be fit from 70 keV up to 2800 keV, and the rear from 300 keV 
to 16 MeV. The huge number of excess counts below 300 keV in the rear detectors is 
understood: The geometrical constellation of the GRB, RHESSI, and Earth was such that 
the GRB photons came from the front direction, where the effective area is relatively small. 
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whereas the Earth was behind RHESSI so that the backscattered photons could easily reach 
the rear segments. 

The only function that fits the front and rear spectra over the entire energy range from 
70 keV up to 16 MeV is the CBM. Fitting front and rear spectra simultaneously with CBM 
and all parameters free, yields (90% CL errors): 



T 


= 678 til keV 




a 


- 60 +°-°^ 

— U.OU ^0.08 


(13) 




9 1 9 +0.08 




b 


= 0.108 tlfi 





with = 187.5 for udoF = 185. These values are used for the black line histogram and the 
residuals in Figs. M and [71 

The front spectrum alone is well described by a Band function. Its is even marginally 
smaller than that of the CBM model fit, see Table [2l The Band function also fits the rear 
spectrum up to 4.5 MeV with = 72.8 (62 DoF), but not at higher energies. The front and 
rear parameters (up to 4.5 MeV) agree, see Table HI Evaluating the full parameter space 
simultaneously for front and rear yields (90% CL errors): 

Ep =71lli^keV 

a = 0.692 ^^39 (14) 
p =3.19 ±0.08 

with = 176.5 for udoF = 174. These values are used for the gr e y line histogram in Figs. [6] 



and [71 They agree with the preliminary results by IWigger et al.l (120071 ). Above 4.5 MeV, a 
PL with 7 = 2.23 ± 0.21 fits the data (x^ = 12.5 for 10 DoF). 



Discussion As can be learned from Table [2l the high energy part can not be fit by Band, 
BPL or CPL, and the low energy part of the spectrum can not be fit by BBPL nor by 
BBmPL. The only function that fits over the whole RHESSI energy range is CBM. 

The CBM function has one parameter more than the Band function. An F-test indicates 
that the chance probability of producing such an improvement in with the additional 
parameter is only 4.0 x 10~^. The spectral hardening at 4.5 MeV is significant. 

Because this GRB has so many counts at high energies, we used a simulation with 
jsim = 1-75 for the results cited above. A power law index of 1.75 results in relatively more 
counts at high energies than the usual power law index (=2). We also used simulations with 
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7sjm = 1.5 and "fsim = 2.0. The results were almost identical, especially for the high energy 
parameters (3 and b of the CBM fit. 

The high energy photon index (3 of the CBM function agrees perfectly with the theo- 
retical expected value (~ 2.1). The low energy photon index a ~ 0.6 on the other hand is 
slightly smaller than expected from theory 1.0). 

Peak resolved analysis The time structure of GRB 021206 is rather intricate. Four 
periods of emission can be distinguished, see Fig. [H each of them probably consisting of 
several overlaying sub-peaks. Luckily, these time periods match quite well our minimum 
time resolution of one third of a rotation period for fitting with a rotation averaged response 
function (see §2.3p . 

The fitted parameters are listed in Table M for the four time intervals marked in the 
figure, as well as the additional tail interval. The tail interval lasts one full rotation, starting 
at the end of the P4 interval. The fluences of the two components in the CBM function 
(eq. [To]) are listed separately {Fcpl for the CPL component and FmPL for the modified PL 
component). The mPL index f] was kept frozen at 2.1. 

The energy Tp increases from the first to the second time interval and then decreases. 
Also FcpL and F^pl increase from the first to the second interval, and then decrease. 
However, the modified PL component seems to decay more slowly than the CPL component. 
The tail is dominated by the mPL component. 



4.4. GRB 030329 



GRB 030329 is famous f or the supernova 2003dh detected in its after glow (iHjorth et al. 



20031 : iMatheson et all l2003l : IChornock et all l2003l : IZaritskv et all l2003f ). The authors of 



the CBM model used the lightcurve of this GRB and its early afterglow to predict the 
later afterglow and the appearance of a supernova fse e iDado. Par &: de Ruiulall2003bl ). A 
supernova an d the l ate afterglow was also predicted by IZeh et al. N2003h . and is discussed in 
Ferrero et al.l (120061 ). 



In the lightcurve of GRB 030329, two peaks are clearly separated (Fig. [H bottom plot). 
We analyse them separately. 



Fit The spectrum of the first peak is shown in Fig. [HI Three bins around a background line 
from 64 keV to 70 keV were not included in the fit. A very good fit (see Table [2]) is a CPL 
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with: Epeak = 158.7 ± 5.0 keV and a = 1.662 ± 0.032. Not surprisingly, the CBM and Band 
function, having more parameters, but being closely related to a CPL, fit only marginally 
better. 

The spectrum of the second peak is shown in Fig. [91 The spectrum has some wiggles 
below 160 keV, which account for the relatively high x^- Many models give an acceptable 
fit, only BBPL does not fit. A good fit (see Table [2]) is a CPL with Ep^ak = 78 ± 13 keV 
and a = 1.83 ± 0.05. Also a good fit is a broken power law (BPL) with Ef, = 175 ± 15 keV, 
a = 1.985 ± 0.032, and (3 = 2.68 ± 0.08. Also Band function and CBM model fit the data 
well, see Fig. |9l 



Discussion The prompt emis sion was detected by HETE. Its spectrura is pu blished by 
Vanderspek et al.l (120041 ) and by lSakamoto et al.l (120051 ). IVanderspek et al.l (120041 ) do a time 
resolved analysis. For the entire burst they find ( 9 0% C L errors): Ep^ak = 70.2 ± 2.3 keV, 
a = 1.32 ± 0.02, p = 2.44 ± 0.08. ISakamoto et"al] J2005h find for the entire burst (90% CL 
errors): Epeak = 68 ± 2 keV, a = 1.26 ± 0.02, /3 = 2.28 ± 0.06. 

The RHESSI parameters, for both peaks, are all significantly higher (see Table H]). 
However, the high energy photon indices can not be compared directly, because the break 
energy (above which (3 is determined, see eq. [7]) for HETE is 116 keV, whereas for RHESSI 
it is > 400 keV, i.e. above the HETE energy range. Fitting the RHESSI data (for the entire 
duration of the burst) from 135 keV to 500 keV only, where the RHESSI response is good, 
we find (3 = 2.441 ±0.032, in excellent agreement with HETE. Fitting the RHESSI data from 
400 keV to 2000 keV, i.e. above the HETE range, we find /3 = 3.11 ± 0.25. The spectrum 
seems to soften above ~ 350 keV. 

The high RHESSI value for a (almost 2.0) for the second lightcurve peak indicates that 
the spectral peak (in the E-Fe representation) is broad. Since RHESSI's sensitivity drops 
below ^ 80 keV, and this is a GRB with Ep^ak in the order of 100 keV, it is likely that 
RHESSI's a describes rather the broadness of the peak than the low energy photon index. 
This opinion is supported by the fit result of the BPL fit. The low energy photon index 
a ~ 2.0 shows that the EEe spectrum is fiat from 34 keV to 175 keV. 



It should also be mentioned that the of the HETE fit, as cited by ISakamoto et al. 



(120051 ). is very bad. Also the RHESSI of the spectral fits for the second peak are rather 
high. A joint fit of RHESSI and HETE data might reveal interesting features. 
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4.5. GRB 030406 

Fit The lightcurve is sliown in Fig. [2] (top) and the spectrum is shown in Fig. [TOl The best 
fit is the CBM function with the parameters a and (5 frozen at the theoretical values. The 
next best fit is the Band function. 



discussion The spectrum of this burst was also studied by iMarcinkowski et al.l (120061 ) 
using data from the INTEGRAL satellite. For the spectral an alysis, they used combined 
ISGRI and IBIS Compton mode data. The time interval used by IMarcinkowski et al.l (120061 ) 
differs from ours. Using a similar time interval as their 'peak' time interval, we find: 

• Fit of a BPL from 24-2400 keV (comparable to the energy range in the analysis by 
Marcinkowski et al.l (120061 )): fits well [y^ = 70.3 for tidof = 60), and the parameters are: 
Eb = 479 ± 80 keV, a = 1.08 ± 0.05, P = 1.96 ± 0.02. 

• Fit of a BPL from 24 keV to 16 MeV: does not fit well (x^ = 86.0 for tidoF = 71). 

• Fit of CBM from 24 keV to 16 MeV with a = 1.0 and (3 = 2.1 frozen: fits well 
(x"^ = 77.9 for udoF = 72), and the parameters are: T = 1220±110 keV and b = 0.055±0.079, 

• Fit of Band function: fits well (x^ = 77.33 for udoF = 71), and the parameters are: 
Ep^k = 1180 ± 120 keV, a = 0.96 ± 0.06, /? = 3.02 ± 0.55. 



For the high energy part, the parameters found by IMarcinkowski et al.l (120061 ) and by 
us agree. But we can not confirm their hard low energy photon index a < 0.0. We even 
dare to say that we trust our low energy photon index better, because for this GRB incom- 
ing direction, the RHESSI response function is well understood, whereas the INTEGRAL 
response function of this burst might suffer from same systematic effects that we described 
for RHESSI in g221 



4.6. GRB 030519B 

Fit The lightcurve is shown in Fig. [2] and the spectrum in Fig. [TTl The best fit is the Band 
function, followed by CBM. 



Discussion In the HETE GRB catalog by ISakamoto et al.l (120051) one finds (90% CL 



errors): Ep = 138 t\l keV, a = 0.8 ± 0.1, /3 = 1.7 ± 0.2. Since f3 < 2.0, the energy Ep is not 
the peak energy, but only a variable related to the parameter Eq = Ep/{2 — a). Indeed, the 
RHESSI peak energy for this GRB is > 400 keV. But the HETE parameters do not fit the 
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RHESSI spectrum from 70 keV to 350 keV [x^ = 141 for 53 energy bins). 

4.7. GRB 031027 

The lightcurve is shown in Fig. [2] and the spectrum in Fig. [121 The best fit is a cut off 
power law with Ep^ak = 336 ± 9 keV and a = 0.940 ± 0.05. This function is shown in Fig. [T^l 
Band function and CBM have difficuhies to converge. Since the CPL fits so well, we expect 
/? = oo for the Band function and 6 = for the CBM function. 

4.8. GRB 031111 

Fit The lightcurve is shown in Fig. [2] (bottom) and the spectrum is shown in Fig. [T31 Band 
function, CBM and BBmPL are good fits, the best being Band function. 

Discussion A preliminary CPL fit to the HETE data is published on a web pag^ as 
Eq = 600.5 keV and a = 0.8366 with a good x^- These values describe the RHESSI 
spectrum well from 80 keV to 350 keV, but not at higher energies. HETE's energy range 
ends at 400 keV, thus we believe that our values for Eq and a are better. 

5. GENERAL DISCUSSION AND CONCLUSION 

5.1. The spectral functions 

What is an acceptable x^? In the limit of many degrees of freedom {tidoF > 30), 
is normal distributed with an expectation value of udof — 0.5 and a variance of a^2 = 
y/2nDoF — 1- A fit is acceptable if is close to its expectation value and if the residuals 
scatter around zero over the whole fit range, i.e. if the fit "looks good". 

From Table [2] we conclude that the CBM gives acceptable for all GRBs studied. And 
they also look good, as can be seen in Figs. [3] to [131 Except for GRB 021206, rear (Fig. [7]), 
the same can be said for the Band function. 

In many cases, a cut off power law (CPL) fits the spectrum up to high energies, e.g. 



3 



http: //space -mit . edu/HETE/Bursts/GRBOSllllA/] 
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GRB 021008, GRB 030329 or GRB 031027. In these cases, Band and CBM improve the 
goodness-of-fit shghtly, but all three spectral shapes fit the data acceptably. 

A broken power law fits sometimes, but usually not well. 

BBPL and BBmPL do not fit in general, BBPL worse than BBmPL. However, it should 
be mentioned that a blackbody com ponent is expected — if at all — only at the beginning 



of a GRB (see e.g. iRyde et al.l ((20061) and references therein), whereas we fitted the entire 



duration of the bursts. When using BBPL, we often find that the PL component fits either 
at high energies or at low energies. This is also discussed by iGhirlanda et al.l (120071 ) who 
studied six BATSE GRBs in detail, where low energy data from the WFC instrument (on 
board BeppoSAX) are available. They find that the WFC data fit the Band function or 
GPL extrapolation, but not the BBPL extrapolation to low energies. Arguing that the PL 
contribution is too simple, they try to fit a blackbody spectrum plus GPL. We suggest to use 
our BBmPL function instead. Its modified PL component describes a spectral break from 
dN/dE oc at low energies to dN/dE oc E"^ at high energies. 



5.2. CBM function 

The present work is, to our knowledge, the first systematic attempt to fit the GBM 
function to prompt GRB spectra. The two terms in eq. [10] have a simple meaning. According 
to the cannonball theory, all GRBs are associated with a supernova. The ambient light is 
Gompton up-scattered by the cannonball's electrons, producing the prompt GRB emission. 
Some electrons are simply comoving with the cannonball, giving rise to the GPL term in 
eq. Uni Since the photon spectrum of the ambient light can be described by a thin thermal 
bremsstrahlung spectrum, a is expected to be ~ 1. The second term (mPL) is caused by 
a small fraction of electror is accelerated to a p ower law distribution, resulting in a ph oton 



index oi (5 ^ 2.1. See e.g. lOado fc Dad (120051 ) . §3.8. or IPado. Par fc de Ruiulal (120071 ). §2 
and 4.1 for a summary. 

In our study, the observed values for a are all approximately 1, as predicted by the 
GBM. Because of the low count statistics at high energies, we could not always fit /?. We 
then fixed it to its theoretical value of 2.1 in order to make the fit converge and to obtain 
a value for the parameter h. In the cases where we could fit /3, we found values close to 2.1 
(Table ED . 

For the factor h of the modified PL component in the GBM function we typically found 
values of the order 0.1 . An exception is GRB 031111, where h is of the order 1.0, but with 
a large error (0.4). 
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Our values for and b are similar to the on es found by iDado fc Pari (120051 ) (fit of 
GRB 941017) and bv lPado. Par fc de Ruiulal tooi l (fit of X-ray flashes XRF 971019, XRF 
980128, and XRF 990520 using BeppoSAX/WFC and CGRO/BATSE data). The authors 
of the CBM hypothesize that XRFs are simply GRBs viewed further off the jet axis. 

The different time development of the GPL- and the mPL-fluences, as reported in Ta- 
ble [6l possibly point to a different time dependence of the two underlying electron distribu- 
tions within a cannonball. 



5.3. Fitting CBM function versus CPL and Band function 

Both the Band function and the GBM function are extensions of a GPL, the Band 
function with one additional parameter, the GBM function with two. For cases where a 
GPL fits the data well, also a Band function with /3 = oo or a GBM function with 6 = 
(and (3 = 2.1 or any other value) fits. This is the case for GRB 031027. 

Whether additional parameters are necessary in a fit, can be tested with the F-test. For 
GRB 030329, the extra parameters are barely needed. For GRB 020715, 021008, 021206, 
030406, 030519B, 031111 additional parameters are required at a confidence level of at least 
90%. 

Goncerning the question of whether the high energy power law parameter f3 in the GBM 
should be treated as free parameter, the answer is 'yes' from a theoretical point of view, but in 
practice, see Table [21 the improvements in are marginal or small for all bursts we studied. 
Our practice to freeze /? at its theoretically predicted value in cases of bad convergence seems 
to be acceptable. 

It is more difficult to compare the goodness of fit using the Band function compared 
to using the GBM function. The two functions are not independent, because they both are 
dominated by a GPL up to the peak energy and higher. In most cases of our study, the 
two functions fit the observed spectrum equally well with a slight preference for the Band 
function. At high energies however (typically above several times the peak energy) the two 
functions are different, the spectral hardening being a unique feature of the GBM function. 
There is only one case, namely GRB 021206, where this hardening is observed. For the rear 
data going up to high energies, a Band function fit gives x^/dof = 133.3/74 (see Table [2]). 
This is not acceptable at < 0.01% probability of being accidentally so high. The GBM fit 
on the other hand gives x^/dof = 82.7/73, which is fully acceptable at a 20% level. 

We would like to stress again that, while the GBM gives acceptable fits for all cases, the 
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Band function fails in one case. This seems enough to us to give some credit to the CBM. 



But it is, of course, no proof that the CBM is the only theory capable of describing the 
spectrum of GRB 021206. For example, a Band function plus a PL with 7 ^ 1.5 would also 
fit. But there is no theory to predict such a shape. To our knowledge, CBM is to date the 
only existing GRB model that explains the prompt GRB spectra from first principles. 

At this place we also would like to note the the mean a- value found for the BATSE 



catalogue is 1 (see iKaneko et al.l l2006l )). We cite from their summary: "M^e confirmed, 
using a much large r sample, that the most common value for the low-energy index is ~ — iQ 



( iPreece et al\\200(A : \Ghirlanda et al\\200B) . The overall distribution of this parameter shows 
no clustering or distinct features at the va l ues expected fr om various emission m odels, such 
as — 2/3 for synchrotron nKatz et al\\l99A:\Tavani\\l99q). for j itter radiation ^Medvedei 
200(a) . or —3/2 for cooling synchrotron ((Ghisellini & Celotti\l99(i) ." They do not mention 
the CBM which would explain a ^ 1. 

Note that the f3 values of the CBM are systematically lower than the (3 values of the 
Band function, compare Tables [3] and |H From Band function fits to BATS E GRBs, it is 



known that f3 is clustered around 2.3, with a long tail towards higher values, see lKaneko et al. 
( I2OO6I ). For CBM we would expect (3 to cluster at slightly lower values. 



For criticism of the CBM, see e.g. lHillad (120061 ) . but see also the answer by lDar fc de Rujula 



(120061). 



5.4. The spectral hardening 

The difference of a CBM spectrum and Band function is the hardening at high energies. 
This becomes visible — for the GRBs studied here — in the few MeV region, but it depends on 
the peak energy and the factor b. For a = 1.0 and b = 0.10 the hardening typically appears 
at several times the peak energy and the second term dominates at 1 times the peak energy. 



For the spectral fit of XRFs done by iDado. Par fc de Rujulal (l2004j ) , the spectral hardening 



is expected in the few hundred keV region, just where the number of photons detected runs 
low. Most of our GRBs also suffer from this lack of statistics at high energies, preventing 
the detection of a hardening. 

A spectral coverage of two decades and good detection efficiency at high energies is 
necessary to experimentally observe the full shape of the CBM function. In the case of 



^ this corresponds to +1 in our notation 
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GRB 021206 we were able to detect this hardening, thanks to RHESSI's broad energy range 
(30 keV to 15 MeV), and because this is one of the brightest GRBs ever observed. 

There is a GRB observed by SMM from 20 keV up to 100 MeV, namely GRB 840805. 
As reported by lShare et al.l (119861 ). the spectrum of this burst shows emission up 100 MeV. 
In order to fit the spectrum, "a classical thermal synchrotron function plus a power law" 
was used. The power law component was required to fit the data above about 6 MeV. This 
is a hint of a spectral hardening around 6 MeV, and we suppose that the spectrum of this 
GRB can be fit by a CBM function. 



The spectral hardening observed in GRB 941017 (iGonzalez et al.l l2003l ) seems to be 
dif ferent. The p h oton index of GRB 941017 above a few MeV is ~ 1.0. This case is discussed 
by iDado Sz Pari (120051 ) as a possible additional feature in the CBM spectrum. 



5.5. Outlook 



In order to find more GRB spectra that show the hardening characteristic for the CBM 
function, strong GRBs have to be observed over a broad enough energy range. With the 
forthcoming GLAST mission, we expect that more such spectra will be observed. But also 
joint analyses with more than one instrument could reveal this hardening. We therefore 
suggest: 

• to search for CBM spectrum candidates among joint Swift/RHESSI GRBs and XRFs, 
and joint Swift/Konus GRBs. 

• to reanal y se so me BATSE bursts. Looking at the BATSE spectra published by 
Ghirlanda et al.l (120071 ). we suppose that the CBM can possibly improve the fits of GRB 
98 0329, GRB 990123 and GRB 990510. The same can be said for GR B 911031 as publishe d 
bv iRvde et aP J2006h . And GRB 000429, as published in Fig. 19 of IXaneko et~all J2006h . 
looks like a promising candidate as well. 



to search in KONUS data for suitable GRBs. 



• to add the CBM function to XSPEC in order to make it more accessible to the 
astronomical community. 
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6. SUMMARY 

We have presented the time integrated spectra of 8 bright GRBs observed by RHESSI 
in the years 2002 and 2003. 

The spectrum of GRB 021206 shows a hardening above 4 McV. From 70 keV to 4.5 MeV, 
the spectrum can be well fitted by a Band function - but not above that. The cannonball 
model successfully describes the entire spectrum up to 16 MeV, the upper limit of RHESSI's 
energy range. For the spectra of the seven other GRBs analysed, we found that they can be 
fitted by the CBM as well as by the Band function. 

We therefore suggest that the cannonball model should be considered for fitting GRB 
spectra. 

We thank K. Hurley, A. Kann, S. McGlynn, J. Ripa and L. Hanlon for helpful discussion 
and comments. 
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Table 1. GRB analysis time intervals 



GRB 


to 


^tburst 


AtBGl 




e 


ref. 




(UT) 


(s) 


(s) 


(s) 


(degrees) 




020715 


19:20:56.0 


[11.53,15.55] 


[ -80.46,0.0] 


[ 48.28,168.97] 


72.4 


1 


021008 


07:00:45.0 


[17.21,21.29] 


[ -73.37,0.0] 


[ 36.68, 48.91] 


50.1 


2 


021206 


22:49:11.7 


[ 2.73, 8.19] 


[ -53.26,0.0] 


[ 20.49,102.43] 


18.0 


3 


030329 PI 


11:37:10.0 


[16.56,24.84] 


[ -70.39,0.0] 


[ 70.39,140.78] 


144.1 


4 


030329 P2 




[28.98,34.50] 










030406 


22:41:30.0 


[85.68,89.83] 


[-140.96,0.0] 


[140.96,281.93] 


96.1 


5 


030519B 


14:04:53.0 


[ 0.46,11.47] 


[ -61.94,0.0] 


[ 28.90, 90.84] 


165.5 


6 


031027 


17:07:06.0 


[29.71,45.92] 


[-137.77,0.0] 


[ 68.88,206.65] 


101.5 


7 


031111 


16:45:12.0 


[ 2.27, 6.35] 


[-122.51,0.0] 


[ 12.25,134.76] 


155.6 


8 



Note. — to'- reference time; Atburst'- time interval for spectral analysis; AtsGi'- back- 
ground time interval before GRB; AtBG2'- background time interval after GRB; time 
intervals are given relative to tp . 9: angle between G RB direction and RHE SSI axis; Ref 



erences: (1) GCN 1456, 1454 flHurlev et al.ll2002ah . (2) GCN 1629, 1 617 flHurlev et al 



2002bh . (3) GCN 1728. 1727 JHurlev et al.l l2002ch. (4) GCN 1997 dVanderspek et al 



20031). {b^ GCN 2 127 flHurlev et al. l2003af). f6\ GCN 2235. 2237 (jLamb et al 



2003; 



Hurlev et al.ll2003cJ) . f7l GCN 2438 JHurlev et al.l[2003dh . (8) HETE trigger 2924, GCN 



2443 (IHurlev et al 



2003d). 
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Table 2. Chi-square of spectral fits 



T> ID 


A LP 

^-C/front 


A JT* 


n 




Band 


CBiVl 






BBFL 


BBmFL 




[KeV ) 


[Ke\ ) 




"^par — 3 


4 


4 


r 




4 


4 


4 


020715 




[ 30,15660] 


117 


113.9 


106.3 


110.8 


110.7 


129.2 


270.0 


157.1 


021008 


[300,2800] 




38 


35.5 


34.8 


35.1 


35.0 


33.0 


32.9 


33.9 


021008 




[300,15660] 


50 


43.8 


39.2 


39.9 


39.2 


52.9 


97.4 


60.2 


021008 


[300,2800] 


[300,15660] 


88 




79.5 


80.8 










021206 


[ 70,2800] 




112 


130.6 


103.6 


104.1 


103.9 


155.5 


315.7 


191.8 


021206 




[300,16000] 


78 


338.1 


133.3 


82.7 


82.7 


132.9 


94.1 


110.5 


021206 


[ 70,2800] 


[300,16000] 


190 








187.5 








021206 




[300, 4500] 


66 




72.8 


69.8 










021206 


[ 70,2800] 


[300, 4500] 


178 




176.5 












030329 PI 




[ 34,10000] 


94 


86.5 


84.3 


84.8 


84.7 


89.3 


137.8 


89.1 


030329 P2 




[ 34, 7000] 


87 


104.5 


103.3 


103.2 


103.1 


98.9 


102.0 


98.4 


030406 




[ 24,15000] 


75 


79.1 


75.6 


75.0 


75.0 


88.6 


151.2 


88.7 


030519B 




[ 70,15000] 


79 


102.2 


86.3 


91.1 


89.0 


99.8 


189.8 


109.8 


031027 




[ 60, 6000] 


63 


63.6 


n.c. 


n.c. 


n.c. 


93.0 


138.9 


99.6 


031111 




[ 38,15000] 


117 


182.4 


128.3 


133.2 


130.4 


140.5 


266.0 


133.0 



Note. — obtained by fitting different spectral models to the data. AEfj.Q^^/j.^g^j.: energy interval used 
to fit front /rear detector data; n: number of energy bins; Hpar- number of free fit parameters; CPL: cut 
off power law eq. O Band: Band function eq. [71 CBM: cannonball model eq. \W\ BPL: broken power law 
eq. El BBPL: blackbody plus power law eq. El BBmPL: blackbody plus modified power law eq. [H n.c: fit 
did not converge. In the case of CBM with 4 parameters, /? was fixed to its theoretically expected value of 
2.1 . For each fit, the degree of freedom is n£)oF = n — Upar- 
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Table 3. Fit results for CBM 



GRB 


/\ rO-P^^vrJ- 

i — 'iront 


A 

^-^-'-^rear 




Oi 





h 
















020715 




[ 30,15660] 


532±20 


0.741 ±0.077 


2.20±0.14 


0.067±0.040 


021008 


[300,2800] 




628±71 


1.31 ±0.28 


2.1 


0.052±0.076 


021008 




[300,15660] 


672±68 


1.487 ±0.062 


2.77±0.55 


0.085±0.092 


021008 


[300,2800] 


[300,15660] 


641±32 


1.523 ±0.055 


2.1 


0.020±0.008 


021206 


[ 70,2800] 




672±20 


0.66 ±0.21 


1.92±0.67 


0.063±0.142 


021206 




[300,16000] 


672±24 


0.67 ±0.19 


2.12±0.13 


0.102±0.048 


021206 


[ 70,2800] 


[300,16000] 


678± 6 


0.60 ±0.06 


2.10±0.08 


0.103±0.028 


021206 




[300, 4500] 


670±23 


0.71 ±0.15 


2.1 


0.091±0.012 


021206 


[ 70,2800] 


[300, 4500] 










030329 PI 




[ 34,10000] 


147±10 


1.614 ±0.036 


2.1 


0.033±0.029 


030329 P2 




[ 34, 7000] 


69±15 


1.841 ±0.049 


2.1 


0.048±0.055 


030406 




[ 24,15000] 


626±83 


0.966 ±0.089 


2.1 


0.18 ±0.12 


030519B 




[ 70,15000] 


396±12 


0.949 ±0.073 


2.388±0.097 


0.135±0.048 


031027 




[ 60, 6000] 


340±17 


0.950 ±0.055 


2.1 


-0.010±0.025 


031111 




[ 38,15000] 


690±45 


0.68 ±0.27 


2.241±0.023 


1.09 ±0.36 



Note. — Ai^front/rear: energy interval used to fit front/rear detector data; Tp, a, /3, h: param- 
eters as defined in the text below eq. [TDl errors are symmetric la errors; where no error is 
given, the parameter was frozen at that value. 
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Table 4. Fit results for Band function 



GRB 


/\ rO-p^^v^j- 
'■ — 'iront 




V 


OL 
















020715 




[ 30,15660] 


531 ± 24 


0.776 ± 0.044 


3.14 ±0.25 


021008 


[300,2800] 




670 ± 58 


1.36 ±0.23 


3.41 ±0.51 


021008 




[300,15660] 


678 ± 42 


1.526 ±0.067 


3.86 ±0.25 


021008 


[300,2800] 


[300,15660] 


677 ± 33 


1.493 ±0.056 


3.73 ±0.18 


021206 


[ 70,2800] 




713 ± 17 


0.694 ±0.031 


3.20 ±0.13 


021206 




[300,16000] 








021206 


[ 70,2800] 


[300,16000] 








021206 




[300, 4500] 


709 ± 18 


0.72 ±0.20 


3.186 ±0.063 


021206 


[ 70,2800] 


[300, 4500] 


711 ±7 


0.692 ±0.020 


3.19 ±0.04 


030329 PI 




[ 34,10000] 


157.2 ±5.2 


1.608 ±0.038 


3.48 ±0.53 


030329 P2 




[ 34, 7000] 


85± 11 


1.781 ±0.065 


3.04 ±0.30 


030406 




[ 24,15000] 


674 ± 70 


0.979 ± 0.064 


2.61 ±0.27 


030519B 




[ 70,15000] 


417 ± 13 


1.048 ± 0.042 


3.11 ±0.18 


031027 




[ 60, 6000] 


338 ± 15 


0.940 ± 0.079 




031111 




[ 38,15000] 


844 ± 59 


1.102 ±0.036 


2.364 ±0.068 



Note. — Ai^front/rear: energy interval used to fit front/rear detector data; Ep, a, (3: 
parameters as defined in the text below eq. [71 where no (3 is given, a CPL (eq. 
was fitted; errors are symmetric la errors. 



-27- 



Table 5. Fluences in 10 ^erg cm ^ 



GRB 




-^Band 


-Frhessi 


-Fhete 


HETE 


-^Ulysses 


Ulysses 










(10~^ erg 








020715 


4.37 


4.41 


3.94 


1.93 




0.43 


0.30 


021008 front 


26.85 


26.85 


41.24 


22.92^ 




8.67*^ 




021008 rear 


35.54 


35.64 


48.85 


27.15^ 




10.27^^ 


8.5 


021206 front 


52.15 


52.44 


55.67 


19.97 




3.81 




021206 rear 


58.74 


53.45 


70.55 


25.29 




4.82 


16 


030329 PI 


6.51 


6.47 


4.42 


5.20 




2.57 




030329 P2 


3.58 


3.58 


2.26 


2.95 




1.69 




030329 total 






7.35 


9.46 


10.76 ±0.14 


4.93 




030406 


4.81 


4.81 


4.26 


1.62 




0.40 


1.3 


030519B 


10.27 


10.36 


9.56 


6.07 


6.10 ±0.1 


1.78 




031027 


5.37 


5.45 


4.81 


3.97 




1.17 


1.4 


031111 


7.40 


7.40 


6.59 


2.10 


1.714 


0.56 


0.21 



Note. — Fcb: fluence from CBM fit (Table[3]); F-Band- fluence from Band function 
fit (Table HD; 

-^RHESsi^ fluence in [100,10000] keV (RHESSI range); -Frete: fluence 
in [30,400] keV (HETE range); Fyiysses: fluence in [25,100] keV (Ulysses range); 
HETE: HETE fluences from references cited in ^ Ulysses: Ulysses fluences from 
references cited in Table [H 

^Our fits overestimate the real counts. More realistic is -Fuiysses= 4.8 x 
10~^ erg cm~^. 
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Table 6. Peak resolved analysis of GRB 021206 



At 


duration 


T 


a 


/5 


h 


FcPL 


FmPL 




(s) 


(keV) 








(10-5 erj 


I cm-^) 


PI 


1.366 


661 ± 18 


0.77 ±0.06 


2.1 


0.059 ±0.010 


11.7 


3.0 


P2 


1.366 


732 ± 12 


0.42 ±0.05 


2.1 


0.115 ±0.010 


14.5 


12.9 


P3 


1.366 


684 ± 14 


0.63 ±0.05 


2.1 


0.085 ±0.012 


12.7 


6.3 


P4 


1.366 


530 ± 20 


0.80 ±0.08 


2.1 


0.113 ±0.018 


6.9 


3.8 


tail 


4.097 


160 ± 60 


1.0 


2.1 


9 c: +00 
^■'J -1.5 


0.2 


3.1 



Note. — At: time period, cf. Fig. [H Tp, a, [3 and h: CBM parameters; Fcpl: 
fluence of the CPL-component in the CBM function (eq. [TD]), -FmPL: fluence of 
the modified PL component of eq. [TOl fiuences are for the range [100,10000] keV. 
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Fig. 1. — lightcurves for the energy band 20 keV to 3 MeV; black: rear detectors, grey: front 
detectors; broken vertical lines: time intervals used for spectral analysis. 
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Fig. 2. — lightcurves for the energy band 20 keV to 3 MeV; black: rear detectors, grey: front 
detectors; broken vertical lines: time intervals used for spectral analysis. 
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Fig. 3. — rear spectrum of GRB 020715; error bars: photon counts after background subtrac- 
tion; black histogram: CBM fit (eq. [TUl) : grey histogram: Band function fit (eq. [7]); vertical 
broken lines: energy range used for fitting; bottom: residuals of CBM (black) and Band 
(grey) fit. 
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Fig. 4. — front spectrum of GRB 021008; explanations see caption of Fig. [3l the same set of 
parameters (eqs. [11] and [T2|) is used for this plot and Fig. [51 
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Fig. 5. — rear spectrum of GRB 021008; explanations see caption of Fig. [3l the same set 
parameters (eqs. [11] and [T2|) is used for this plot and Fig. |H 
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Fig. 6. — spectrum of GRB 021206, front detectors; explanations see caption of Fig. [3l the 
same set of parameters (eqs. [13] and [T^ is used for this plot and Fig. [71 
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Fig. 7. — top: spectrum of GRB 021206, rear detectors; symbols as in Fig. [3l bottom: 
residuals of CBM fit (black) and Band fit (grey). The Band function was fitted only up to 
4.5 MeV. The excess counts below 300 keV are backscatters from Earth, see text. The same 
set of parameters (eqs. [13] and [H]) is used for this plot and Fig. [61 
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Fig. 8. — rear spectrum of GRB 030329, first peak; explanations see caption of Fig. [31 
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Fig. 9. — rear spectrum of GRB 030329, second peak; explanations see caption of Fig. [31 
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Fig. 10. — rear spectrum of GRB 030406; explanations see caption of Fig. [31 
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Fig. 11. — rear spectrum of GRB 030519B; explanations see caption of Fig. [31 
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Fig. 12. — rear spectrum of GRB 031027; histogram and residuals: fit of a cut off power law 
(CPL); a CPL is equivalent to a Band function with /? = oo or CBM with 6 = 0. 



-41 - 



> 



> 
CD 



c 

O 

o 



o 

(D 




100 1000 
energy [keV] 



10000 



Fig. 13. — rear spectrum of GRB 031111; explanations see caption of Fig. [31 



